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Light hadron production in inclusive pp-scattering at LHC 
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The inclusive production of light mesons in pp-scattering is considered in the framework of reggeon 
phenomenology with supercritical Pomeron. Available low-energy data can be explained with three 
reggeon particles taken into account. With the results obtained rapidity and pseudorapidity distri- 
butions for light-meson production at the LHC energies are predicted. 
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I. INTRODUCTION 



O I It is well known that the growth of total cross sections can be explained in the framework of reggeon theory with 

^^ ] supercritical Pomeron [l| . In this approach the amplitude of an elastic hadron-hadron scattering is written down as a 

Q^ . sum of reggeon diagrams which corresponds to multiple rescattering. The first term of such eikonal expansion is the 

pvi ■ contribution of a "bare" Pomeron: 



P,; A{s,t) = ^^^-cot^js'^^*>^{t), (1) 

QJ ' where the Regge trajectory in the j-plane is parameterized by a linear function a{t) — 1 + A + a't with the intercept 

^ '. a(0) = l + A. 

According to the perturbative QCD calculations [2] the vacuum exchange singularity has rather complex form. It 

r/-\ , can be expressed as a sum of poles in the angular momentum plane with 1 < j < 1 + A, where 
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For as = 0.2 H- 0.25 this allows A to approach rather large value of 0.55 ~ 0.65. 

The information on the parameter A can also be obtained from the analysis of deep inelastic scattering structure 
functions in the region of small Bjorken variable x |3|, where the Regge asymptotic is valid. The power parameter A, 
defined as A = d\nf/d{l/x) - 1, varies from A = 0.23 at Q^ = 16GeV^ to A = 0.5 at Q^ = 400 GeV^. 

The data on total cross sections were analyzed using eikonal model with different values of parameter A. In the 
earlier works [J| total cross sections were fitted using A ^ 0.13. In the paper |5| a double-pole approximation for bare 
Pomeron was used for the first time with A being in the region 0.09 < A < 0.3. In the recent work [a| the triple-pole 
^ ' approximation was used with Ai = 0.058, A2 — 0.0167, A3 = 0.203. 

Other approach to determine the parameters of bare Pomeron is connected with the inclusive production in the 
central reg ion [3-l3]- The argument in support of such method is the AGK cancellation of cut contributions in this 
region [10|. Therefore one deals only with the first term — double bare Pomeron diagram. As a result, the asymptotic 
of inclusive production cross section in the central region has a simple form: 



da 
dy 



(2) 
y=o 



As it was mentioned earlier, only the leading term is taken into account in the expression ([2]), while the contributions 
of secondary trajectories are neglected. However, at low and moderate energies (^/i < 100 GeV) the growth of particle 
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density is caused by contributions of secondary trajectories in Mueller-Kanchelli diagrams [ll[ . These terms are also 
important at higher energies for large values of rapidity y ~ Ina; + In {y/s/m±). The analysis of experimental data on 
da/dy\y=o presented in our previous works [1,|0| shows that power behavior of this inclusive cross section is caused by 
the Pomeron singularity with A « 0.17. In that analysis only low-energy data with y/s < 900 GeV were used. In the 
present work we are going to check whether new measurements of cross-sections of inclusive production of charged 
particles and ii'5-mesons, obtained at LHC at y/s — 0.900, 2.36 and 7 TeV energies agree with predictions based on 
the low-energy data analysis. 

II. DOUBLE-REGGEON EXCHANGE 

According to the generalized optical theorem and factorization of the leading regge singularities the inclusive cross 
section of the reaction a5 — > c + A" in the central region can be written down as 

dV _ 1 d^a _ 1 /_t\"'(0) /_u\"^(°) 

d^Pc IT dni^dy s ^—^ V ^o / \ so / 

where «,j — R, P and F for Reggeon, Pomeron and Froissaron exchange respectively. The intercept parameters for 
these particles are 

Afl = -0.5, Ap = 0.06, Af = 0.17. 

The kinematical variables are defined as t = (pa — Pc) ~ —^/smxe^^, u = (pt, — Pc) ~ —\fsmTe^., s = [pa + Pb) , 
niT = \/ti'v^c +Pt7 So = 1 GeV , and rapidity y is defined as 

Ec^PL 



y = lU' 



rriT 



where E^, pl and pt are energy, longitudinal and transverse momentum of particle c in the laboratory frame. The 
dependence on transverse momentum is contained in the vertex functions fij{mT) which are determined from the 
experimental data. In our work the following parametrization for these vertex functions is used: 

ftj (rriT) = A^jipij (ttt-t) , 

where constants Aij are determined from the fit of low-energy data (-^/s < 900 GeV) on cross-sections of inclusive 
pp ^ c + X reactions in the central region. Transverse momentum dependence of (/^^(tot) -functions can not be 
determined from Regge phenomenology, so, following the work T], a simple exponential parametrization is used in 
our article: 

This functions obey a normalization condition 

00 
ipij {ttit) drrirp = 1. (5) 

Parameter (ipp for Pomeron singularity was determined from the fit of low-energy data and equals 

13 pp = 6GeV-^ 

For A = 0.17 singularity, which dominates at high energies, pT-spectrum of charge particles at LHC energies was 
used to obtain 

Pff = 4.7GeV~^ 
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TABLE I: Numerical values of parameters Aij for c and Ks production in inclusive pp-scattering 

The decrease of this slope parameter with the growth of A can be connected with the decrease of the effective 
radius of the Pomeron singularity, observed, for example, in the description of total and elastic cross sections in the 
three-Pomeron model [y| . 

Let us first consider rapidity distributions at various energies squared s. Integrating the expression ([3]) over the 
transverse mass one obtains 
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dy 
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cosh[(A, -Aj)y]. 



It is clear that this distribution does not depend on slope parameters I3ij. It is important to note, that the in 
presented expression y- and -yi-dependencies are strongly correlated. Thus the rapidity distributions at fixed ^/s are 
as important as the -^/i-dependence at fixed y. 

Experimental data on these distributions and references to original works are given in Q (see references [10-27] in 
that paper). Using this data we obtain the values of parameters A^j, presented in Tab. HI 

In Fig. [T] rapidity distributions of inclusive production of charged particles (upper figure) and Ks-Taeson (lower 
figure) in proton-proton interactions at different energies and energy dependence of da /dy cross section at y = 
are given. Dots in these figures represent experimental data used in the fitting procedure. Theoretical description, 
obtained in the framework of Regge model with the parameters Aij presented above, are shown with lines. 

Experimental data on charged particles distributions at LHC energies are usually reported as distributions over 
pseudorapidity 
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where 9 is scattering angle of the final particle. The relation that connects rapidity and pseudorapidity is 



u = - In 
'^ 2 



pip cosh 1] + pt sinh -ij 



pip cosh T] — Pt sinh 77 



It is clear, that in the limit of zero mass one has rj = y. 

To transfer from rapidity to pseudorapidity distribution the information on the WT-dependence of vertex functions 
(fiij is needed: 



da 
drj 
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E dydrri^ 
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It was earlear mentioned, that Regge theory gives no information about the form of these functions. So the exponential 
parametrization (|4]) which agrees with the experimental data is used. In addition, we assume factorization of the slope 
parameters of these functions, so that the following relation is fulfilled: 



where 



A + /3, 



/3ii = 2.5GeV"\ /3p = 3GeV"\ /3f = 2.35 GeV^^ 
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FIG. 1: Fits of available experimental data by double-Reggeon exchange model: 

(a) — rapidity distributions of pp — ^ c~ + X reaction cross section. Curves and data sets in the figure correspond to energies 
^/s = 4.9 GeV, 11.5 GeV, 13.9 GeV, 21.7 GeV and 53 GeV from bottom to top; (b) — energy dependence of differential cross 
section of the reaction pp ^ c~ + X at y = 0; (c) — rapidity distributions of pp — >■ Ks + X reaction cross section at ^/s — 13.76 
GeV, 63 GeV and 546 GeV; (d) — energy dependence of the differential cross section of the reaction pp — >■ Ks -\- X a.t y = Q 

Now we have enough information to obtain pseudorapidity distributions of the cross sections oi pp ^^ c + X reaction 
at different energies. Experimental data at LHC energies are reported normalized to non-single diffraction cross 
section ctnsd • The following parametrization for aNSD is used in our work: 



o'nsd(s) 



1.76+19.8 



iGeV2 



0.057 



nib. 



In Fig. [5] (a) pseudorapidity distributions at high energies are shown in comparison with experimental data taken 
from the works [13l4l5{ . In these works distributions for both positively and negatively charged particles are given, 
thus they are doubled in comparison with those presented in Fig. [TJ Energy dependence of differential cross section 
of the reaction pp ^ c^ + X &i zero rapidity is presented in Fig. [2] (b) together with both low-energy data taken from 
Q and recent LHC results [l3l - [la |. One can clearly see, that in the framework of Regge phenomenology we get good 
agreement with the higher energy experimental data, although these values were not used in the fitting procedure. 

As ii's-mesons have definite mass, distributions over rapidity y are experimentally measured. Recent data from 
the LHC contains only LHCb results in 2.5 < y < 4 region at ^/s = 0.9 TeV. These points together with the Regge 
model predictions at ^/s — 0.9, 2.36, 7 and 14 TeV are shown in Fig. [3] (a). The agreement within the error bars 
is achieved but experimental points are systematically below the theoretical prediction. These measurements are 
reported in units of cross section as opposed to others, reported in units of multiplicity. Thus the uncertainty from 
the normalization cross section may arise. The rapidity density at y = have been measured up to 1.8 TeV energy 
|16l4l9l |. All these measurements are in pretty good agreement with our theoretical predictions. 
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FIG. 2: (a) — pseudorapidity distributions of pp — ?> c + X reaction cross section. Curves in the figure correspond to energies 
-^s — 0.9, 2.36, 7 and f4 TeV from bottom to top; (b) — energy dependence of differential cross section of the reaction 
pp — >■ c + X at ?7 = 0. Recent values obtained by ALICE and CMS collaborations are marked with stars. 
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FIG. 3: (a) — rapidity distributions of pp — > Ks + X reaction cross section. Curves in the figure correspond to energies 
V^ — 0.9, 2.36, 7 and 14 TeV from bottom to top; (b) — energy dependence of differential cross section of the reaction 
pp — >■ Ks + X at r; = 0. 



III. DISCUSSION 



The energy dependence of total and inclusive cross sections of particle production in central region is usually 
discussed in the framework of Regge phenomenology. In this approach the amplitude of the process is written as 
a sum of reggeonic diagrams: Pomeron and secondary trajectories. According to perturbative QCD the Pomeron 
singularity can be presented as a number of poles in angular momentum plane with 1 < j < 1 + A. 

It was sufficient to consider two Pomeron poles with the intercepts q(0) = 1 and a(0) — 1.17 to describe low-energy 
data on proton-proton scattering, published long before LHC start. This approach gave both energy dependence 
of total cross sections and rapidity distributions at different energy values [9*]. We would like to note that these 
distributions are strongly correlated with each other. The constants of Pomeron interaction with protons and final 
particles can be determined from a simple fit of experimental data. In our article we use low-energy data at ^/s < 900 
GeV energies to determine Pomeron coupling constants and predict rapidity and pseudorapidity distributions of 
light charged and i^g-mesons at LHC. At ^/s — 0.9, 2.36 and 7 TeV energies good agreement with the available 
experimental data is observed. We also give predictions of pseudorapidity distributions of charged particles and 
Ks-Taeson production at ^/s — 7 and 14 TeV. 
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